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ABSTRACT 

We present the results of a study of the structure of a parallel compressible mixing layer in 
a binary mixture of gases. The gases included in this study are hydrogen (H 2 ), helium (He), 
nitrogen (N 2 ), oxygen (0 2 ), neon (JVe) and argon (At). Profiles of the variation of the Lewis 
and Prandtl numbers across the mixing layer for all thirty combinations of gases are given. It is 
shown that the Lewis number can vary by as much as a factor of eight and the Prandtl number by 
a factor of two across the mixing layer. Thus assuming constant values for the Lewis and Prandtl 
numbers of a binary gas mixture in the shear layer, as is done in many theoretical studies, is 
a poor approximation. We also present profiles of the velocity, mass fraction, temperature and 
density for representative binary gas mixtures at zero and supersonic Mach numbers. We show 
that the shape of these profiles is strongly dependent on which gases are in the mixture as well 
as on whether the denser gas is in the fast stream or the slow stream. 
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1 Introduction 


In the last decade the structure and stability of compressible mixing layers has become, 
once again, an area of very active research. There are many experimental studies (Brown 
and Roshko^; Bogdanoff^; Papamoschou and Roshko®; Samimy and Elliot^; Elliot and 
Samimy^; Clemens and Mungal®; Hall, et al.l^; Elliot, et ah®) of the overall structure and the 
turbulence, principally the large scale structures, of these flows. In some of these experiments, 
for example those of Papamoschou and Roshko®, Hall, et al.^ and Ninnemann and Ng®, 
different gases are in the two streams. In other words, the mixing layer exists in a binary gas 
mixture. 

There are also a substantial number of theoretical and numerical studies of the structure of 
the compressible mixing layer as well as the stability of this flow. In almost all cases, it was 
assumed that the fluid properties in the mixing layer were those of a single gas and that variations 
of density and viscosity across the layer were due solely to the variation in the temperature. 
Of course as the mixture fractions of the gases change across the mixing layer, the viscosity, 
diffusivity, specific heat, etc of the mixture will change so as to reflect the changing composition 
of the mixture as well as the change in the temperature. 

Jackson and Grosch^ formulated the stability problem of a compressible mixing layer for a 
gas mixture, assuming that the mean flow field was known. They showed that the dimensionless 
sonic speeds of the two streams as well as the convective Mach number were dependent on the 
composition of the gas mixture. Kennedy and Gatski^^ carried out a theoretical calculation of 
the structure of a supersonic/supersonic mixing layer in the binary gas mixture resulting from 
the mixing of a hydrogen stream with a nitrogen stream. Calculations were carried out for five 
values of the convective Mach number. The Kennedy-Gatski study included the variation of the 
fluid properties across the mixing layer caused by the variation of the gas composition. The 
purpose of this paper is, in part, to correct and, in part, to extend to a wider variety of binary 
gas mixtures the study of Kennedy and Gatski^^. In particular, we will show that certain of 
their results are incorrect and certain of their conclusions must be modified. 

The gases included in this study are hydrogen (# 2 ), helium (He), nitrogen (JV 2 ), oxygen 
(O 2 ), neon (Ne) and argon (Ar). The choice of gases was not arbitrary. Hydrogen is the 
fuel for proposed scram-jet engines. Papamoschou and Roshko® and Hall, et alJ^ used binary 
combinations of He, N 2 and Ar in their experiments. An accurate calculation of the structure of 
flow field in a gas mixture and a meaningful comparison of the results of theory with experiment 
must incorporate not only the variation of density with composition but also the changes with 
composition and temperature of the thermal properties and the transport properties of the 
mixture . 

2 Thermal Properties and Transport Coefficients 

The binary gas mixture is characterized by the mass fractions, { Fj} , j = 1, 2 of the species. 
The mixture properties, heat capacity (Cp), coefficient of viscosity (p), thermal conductivity (k) 
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and binary diffusion coefficient (Z?i 2 ) are functions of the {Fj} and of the temperature T . These 
properties vary continuously through the mixing layer as the temperature and mass fractions 
vary. 

The quest for an understanding of the variation of these thermal and transport properties 
with temperature and mixture composition has been the central topic of the Kinetic Theory of 
Gases and has engendered hundreds, if not thousands, of experiments for well over a century 
and a half (Chapman and Cowling! Hirschfelder, et a]J^; Touloukian, et alJ^b Maitland, 
et al.f^). Indeed, Chapman and Cowling’s book, “The Mathematical Theory of Non-Uniform 
Gases", is subtitled “An Account of the Kinetic Theory of Viscosity, Thermal Conduction and 
Diffusion in Gases”. Despite the vast effort which has been expended on this topic, the theory 
is far from complete and yields some seemingly paradoxical results. 

A satisfactory theory for the variation of transport properties of monatomic gases, both in the 
unmixed state and as a mixture, has been developed (Chapman and Cowling^; Hirschfelder, 
et al.^1; Touloukian, et alJ^; Maitland, et al.^) in terms of the molecular scattering cross 
sections, temperature and composition. For binary mixtures the theory is reasonably straight- 
forward but becomes much more complex for mixtures of more than two gases. There is good 
agreement between theory and experiment for binary gas mixtures. However, there is much less 
experimental conformation of the theory for other than binary mixtures. 

The theory of the transport properties of mixtures of monatomic and polyatomic gases and 
of binary mixtures of polyatomic gases is in a far less satisfactory state. The basic kinetic 
theory for the transport properties of a pure polyatomic gas is a semi- classical one with certain 
quantum mechanical corrections (Touloukian, et al.^; Maitland, et al.!^). It appears that 
there is a great sensitivity of the results to the detailed structure, including the anisotropy, of 
the molecular potentials of polyatomic gases and these are not well known. In addition, some 
of the collisions are inelastic or partially elastic and this adds further difficulties. For mixtures, 
the theory is even more complicated and is, at best, no more accurate than the first order 
Chapman- Cowling approximations and usually less accurate because of other approximations 
which must be made in order to make the theory useful. The theory must be supplemented by 
experimental measurements, sav of the variation of viscosity with temperature, in order to obtain 
the values of intermolecular potentials required to calculate the other transport coefficients. 
These experimental measurements are, in many cases, lacking or inadequate. 

In the absence of a complete theory, we have used a mixture of the results of experimental 
measurements and approximate theory to obtain the thermal properties and transport coeffi- 
cients of the binary gas mixtures. 

2.1 Heat Capacity 

The thermodynamic properties of a large number of substances over a wide range of tem- 
peratures are tabulated by McBride, et alJ^. The heat capacity of the individual gases used in 
this study were taken from this tabulation and stored, as a function of temperature, in tabular 
form within the code. The values for the individual gases, as a function of temperature, were 
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obtained by interpolation. The mixture Cp was calculated from : 


C P = F x C Px + F 2 Cp 2 . 


2.2 Coefficient of Viscosity 

The viscosity of the individual gases was calculated using the first order Kinetic Theory 
formulae as given bv Maitland, et alJ^ (see their Appendix A5.2) and also by Anderson^^. 
The viscosity is an explicit function of temperature and of a collision integral which, in turn, is 
also a function of temperature. Maitland, et alJ^ tabulate values of various collision integrals 
(see their Appendix 2) as a function of temperature. We used the values for the Lennard- Jones 
(12-6) potential given in table A2.2. The calculated viscosity of the individual gases was stored 
within the code in tabular form as a function of temperature. The value of the viscosity at a 
specific temperature was obtained by interpolation. 

The viscosity of a binary mixture of monatomic gases is given by Maitland, et al.^ as a 
function of the mole fractions of the gases, the viscosities of the individual gases and of a so-called 
interaction viscosity arising from collisions between the molecules of the different gases. We 
converted the mole fractions to mass fractions and used the values of the Lennard- Jones (12-6) 
collision integrals in the calculation of the interaction viscosity. Again, the values of interaction 
viscosity were stored within the code in tables as a function of temperature. Interpolation was 
used to obtain the value at a specific temperature. The value of the mixture viscosity at a 
specific temperature was obtained using the formulae of Maitland, et alJ^l with the {F ; } of 
the individual gases, their viscosities and the interaction viscosity. 

It should be noted that the theory used was that for a binary mixture of monatomic gases. 
As pointed out by Maitland, et alJ^ (see page 315 and following) at first order the form of 
the equations for the viscosity of a binary gas mixture of polyatomic gases is the same as for a 
binary mixture of monatomic gases. For this reason, we used the formulae for monatomic gases 
to also calculate the viscosities of polyatomic gases. Comparison of the calculated viscosities for 
the polyatomic gases with tabulated experimental results (Touloukian, et alJ^l) showed good 
agreement; the differences ranged from 0.01% to 0.20%. 

2.3 Thermal Conductivity 

We attempted to use a theory, similar to that used for viscosity, to calculate the thermal 
conductivity. The calculated values were in very poor agreement with tabulated experimental 
data for the polyatomic gases ( H 2 , N 2 and 0 2 ). This lead us to use the tabular experimental 
data (Touloukian, et alJ^) for the thermal conductivity of the individual gases. The mixture 
conductivity was then calculated using the method of Mason and Saxena recommended by 
Touloukian, et alJ^. Because this theory has one form for mixtures of monatomic gases and 
one for mixtures of polyatomic gases, the polyatomic form was used if at least one of the gases 
was polyatomic. 
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2.4 Binary Diffusion Coefficient 

Maitland, et give a formula from first order Kinetic Theory for the binary diffusion 

coefficient, Du of a gas mixture. This is a function of temperature and involves a collision inte- 
gral which is also a function of temperature. In contrast to mixture specific heat, viscosity and 
thermal conductivity, D \ 2 is independent of the mass fractions of the gases. There is relatively 
little experimental data for Du- We compared the computed values with the experimental val- 
ues given by Hirschfelder, et alJ 13 i and by Bzowski, et alJ 19 l The agreement was found to be 
satisfactory; in most cases differences were in the range of 1 % to 2 % with the maximum differ- 
ence being 3 . 3 %. As with the other coefficients, Du was computed at a series of temperatures 
for each combination of gases. These values were stored in tabular form and values at specific 
temperatures were found by interpolation. 


3 Basic Equations 


In what follows in referring to a binary gas, Gas\ — Gas 2 , we use the convention that Gas\ 
has a mass fraction F\ which is zero at +oo and that Gas 2 has a mass fraction F 2 which is zero 
at -oo. That is, the stream at -oo is entirely Gasi and that at +oo is entirely Gas 2 . These 
gases have molecular weights Wi and W 2 , with W = W 2 /W x being their ratio. The velocity 
field has components (U, V), T is the temperature and p is the density. 

The equations were made nondimensional by the values at +00 of the density, poo, tempera- 
ture, T^, speed, Uoo, pressure. = TZpoaT<x,/W 2 ( TZ is the universal gas constant ), viscosity, 
Poo, thermal conductivity, mass diffusion coefficient, £>12,00 , and specific heat, Cp i00 . The 
ratio of specific heats at +oo is and Moo = CW a oo is the Mach number at +00 with the 
speed of sound given by 

a lo = (Too - l)C>oo T’oo- 


The Prandtl and Lewis numbers are defined as usual 


Pr = 


K 


( 1 ) 


and 


Le 


K 

p DuCp 


( 2 ) 


These vary across the shear layer because the transport coefficients and thermal properties vary 
across the layer. 

Using the Howarth-Dorodnitsyn transformation, 


Y = 



V = P V + U 



px 


dy , 


( 3 ) 


and the boundary layer approximation, 
similarity variable, 


with primes denoting differentiation with respect to the 



( 4 ) 
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and with 


u = f(tj), V = (77 f - f)/y/x, (5) 

the nondimensional equations for the steady, laminar flow of compressible binary mixture are. 



ip? /")' + 2 //" = 0, 

(6) 

((£)' or) 

) + 2 Cp / r + (700 - 1) Ml Pfi (/" ) 2 + (p^) T' J2 CpjF ’ = 0. 

(7) 


(UieH) +2 ^ = 0 ' 

(8) 


1 = F\ + Jp2- 

(9) 


1 = pT (WF l + F 2 ), 

(10) 


C P = C Pl F, + Cp 2 F 2 . 

(11) 

The boundary conditions are 



/' = T = F 2 = 1, Fi = 0 77 — ► + 00, 

(12) 


f — flu , T = (3t , F\ = 1, F 2 = 0 ij — ^ - 00. 

(13) 


The parameters 0u £ [ 0 , 1 ) and fij > 0 are the ratios of the speed and temperature of the gas 
stream at —00 to that at +00. If f 3 j is less than one the gas in the slow stream is colder than the 
gas in the fast stream. If 0t is greater than one the converse is true. The remaining parameter 
is W, the ratio of molecular weights. If W > 1 the gas in the fast stream is heavier than that in 
the slow stream, and if W < 1 the gas in the fast stream is lighter than that in the slow stream. 

The equations given by Kennedy and Gatski^^ differ from those above because of their 
inconsistency in defining the Lewis number. Their definition of the Lewis number (just after 
equation (5) on page 664) is the same as that used here and they used this definition in calculating 
the Lewis number. However, in their equation (5) the Lewis number is the reciprocal of that 
defined. Their equations (18) through (22) perpetuate this inconsistency. In addition, they used 
the kinetic theory for mixtures of monatomic gases (even though both H 2 and JV2 are polyatomic) 
to calculate the transport coefficients. They then used an ad hoc scheme to correct the thermal 
conductivity to account for polyatomic effects. In contrast, we here use a combination of the 
theory for polyatomic gases and tabulated experimental data. 

The system of equations for /, T and F\ is seventh order and there are only six boundary 
conditions. A seventh boundary condition can be found by matching the free stream pressures 
across the layer (Ting^O]) if at least one of the streams is supersonic. However, if both are 
subsonic, Klemp and Acrivos^ showed that this condition is incomplete. The seventh boundary 
condition is equivalent to determining the value of the stream function / at 7? = 0 . We have set 
/( 0 ) = 0 . 

The system of equations and boundary conditions is a two point boundary value problem 
which can be solved by a shooting method. In the limit ?? — — oc the asymptotic forms of 
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(6) - (8) can be integrated exactly. Using the three boundary conditions for /'(-oo ), T(-oo) 
and Fi(-oc) to determine three of the constants of integration, these solutions contain four 
unknown constants. Choosing values for these four constants determines the complete initial 
conditions for (6)-(8) at r/ - -oo. A fourth order double precision Runge-Kutta method was 
used to integrate the system to rf — oo. Iteration on the values of the four constants was carried 
on until the condition at rj = 0 and the boundary conditions at 77 = 00 are satisfied. The system 
was found to be extremely stiff and a nested iteration scheme was used. In this scheme, one 
condition at a time was used and after each was satisfied additional iterations were required to 
resatisfy the previously used conditions. This iterative process was continued until all of the 
boundary conditions w r ere satisfied. 

4 Results 

Kennedy and Gatski^^ only considered the mixture N 2 -H 2 ■ Here we considered the thirty 
different combinations of H 2 , He , Ne, N 2 , 0 2 and Ar for Gas\ - Gas 2 . We did not consider 
systems of the same gas in both streams because this would require the inclusion of self-diffusion 
for consistency and the theory and experimental measurements for this effect are wanting. Table 
1 contains the values of W for each combination. This parameter has a wide range, from 19.813 
for a H 2 — Ar mixture to 0.050 for a Ar - H 2 mixture and there is a corresponding variation 
in the structure of the mixing layer. We first present the variation of the Lewis and Prandtl 
numbers across the shear layer at M 00 = 0. In all cases we took T ^ = 300°/i', fiu = 0.5 and 
3t = 1 . 5 - 

Figures 1, 2 and 3 are plots of the variation of the Lewis number across the mixing layer. 
For most pairs of gases there are relatively large variations in the Lewis number. The largest 
variation occurs for mixtures of He with Ar, 0 2 and N 2 with the ratio of the values at ±00 for 
He — Ar and Ar - He being about eight. Large variations in the value of the Lewis number 
also occur for mixtures of H 2 with Ar, 0 2 and iV 2 . It should be noted that the magnitude of 
the variation is not strictly monotone with molecular weight ratio. The values of W are larger 
for combinations of H 2 with the heavier gases than for He and heavy gases. Yet the variation 
of the Lewis number is smaller for combinations of H 2 with heavy gases than for He with heavy 
gases. A few gas mixtures have nearly constant Lewis number; mixtures of 0 2 , iV 2 and Ar as 
well as the H 2 — He and He - H 2 mixtures. It might seem that a nearly constant Lewis number 
would occur only for mixtures with W close to one. This is not true. The molecular weight 
ratio for mixtures of 0 2 , N 2 and Ar lies in the range of 0.8 to 1.2, and these mixtures have 
relatively constant Lewis number across the layer. However, mixtures of of tf 2 and He have W 
close to 2.0 (or 0.5) and still have small variation of the Lewis number across the layer. Finally, 
it can be seen that the Lewis number profiles are generally not symmetric about the centerline 
of the mixing layer and that interchanging the gases between -hoc and — oc does not yield mirror 
images about the mixing layer centerline. Particularly noticeable cases are mixtures of H 2 or 
He with any of Ar, A r e, N 2 or 0 2 - 

The variation of the Prandtl number across the mixing layer is shown in Figures 4, 5 and 
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6 a. 6 b for all of the gas mixtures. Figure 4 is a plot of the Prandtl number for the same gas 
combinations as those whose Lewis number variation is shown in Figure 1. Figures 5 and 2 have 
a similar relation. However, for the gases whose Lewis number variation is shown in Figure 3, 
a single figure could not adequately display the Prandtl number variation. As a result, Figure 
6 has two parts (a) and (b). There is a substantial variation of the Prandtl number across the 
mixing layer for all combinations of H 2 or He with the heavier gases Ar , Ne, N 2 or 0 2 . As with 
the Lewis number variation, the Prandtl number curves are generally not symmetric about the 
center of the mixing layer, nor does interchanging the gases between +00 and -00 yield mirror 
images of the Prandtl number curves. 

In comparing these results to those of Kennedy and Gatski^ 1 ! there is general disagreement 
between the results for the Prandtl and Lewis numbers given here and those of Kennedy and 
Gatski. The shapes of the Prandtl and Lewis number profiles given here are, approximately, 
“mirror images” about rj = 0 of those given by Kennedy and Gatski. The Prandtl number is 
0.700 in the slow jV 2 stream, but only about 0.655 at g = - 20 . 0 . Its value decreases slowly 
towards the center of the shear layer, reaching a minimum of about 0.38 just short of the middle 
of the shear layer and then rises very rapidly with distance beyond the center of the shear layer, 
reaching a value of 0.711 in the fast H 2 stream. The Lewis number distribution is, on the other 
hand, a monotone function of distance through the shear layer. It rises slowly from a value 
of 0.322 in the jV 2 stream until near the center of the shear layer and then asymptotes over a 
short distance to 2.010 in the H 2 stream. The different shape of the profiles of these parameters 
between our results and those of Kennedy and Gatski is due to the inconsistent definitions used 
by Kennedy and Gatski as discussed in Section 3. 

Finally, the lack of a reflection symmetry of both the Lewis and Prandtl number curves about 
the center of the mixing layer is a result of the non-symmetric variation of the temperature across 
the layer. As stated above, all of the transport coefficients are functions of temperature and so 
show a corresponding non-symmetric variation across the layer. The “tails” of the Lewis number 
and Prandtl number curves are always longest on the side of the mixing layer where a lighter 
gas is diffusing into the higher temperature region as compared to the reverse situation. A good 
example of this is provided by the Lewis number and Prandtl number curves for Ar — H 2 as 
compared to those for H 2 — Ar. 

We next present additional results for the case considered by Kennedy and Gatski^^, N 2 —H 2 
for which W = 0.072. For this case the slow speed stream consists of N 2 and the high speed 
stream is H 2 . Figure 7 shows profiles of the dimensionless velocity, mass fractions, temperature 
and density as a function of the dimensionless distance across the shear layer for M ^ = 0. The 
general shapes of the velocity, mass fraction, temperature and density profiles are similar to 
those obtained by Kennedy and Gatski. This is, in part, because the values of these variables at 
± 00 are set by the boundary conditions and the equations enforce a smooth transition between 
these boundary values. As seen from Figure 7, the velocity profile has a modest “tail” into the 
slow stream and a rapid transition to the fast stream. The mass fraction of H 2 decreases quite 
slowly on the slow speed side of the mixing layer (77 — * - 00 ); there is substantial diffusion of the 
H 2 from the relatively cold fast stream into the the relatively hot slow stream. This is also clear 
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from the Lewis and Prandtl number profiles shown in Figures 3 and 6(a). This effect is also 
apparent in the temperature and density profiles which have very long ‘‘tails" on the N 2 side of 
the layer; thus the temperature and density profiles are shifted towards the N 2 side of the layer 
as compared to the velocity profile which lies almost exactly in the center of the layer. The 
mid-points of the mass fraction profiles are shifted towards the N 2 side of the layer although not 
as much as the density and temperature profiles. Our results are similar to, but not identical to 
those of Kennedy and Gatski. We also find that one of the points made by Kennedy and Gatski 
is still valid; assuming constant values for the Lewis and Prandtl numbers in the shear layer is 
a poor approximation when considering binary gases. 

Next, the corresponding profiles, for M ^ = 0, in the case of H 2 — A 2 {W = 13.897) is shown 
in Figure 8. In this case the slow speed stream is H 2 and the high speed stream is N 2 . Note 
that Kennedy and Gatski^ ^ did not present any results for this case. The mean profiles of 
velocity, mass fraction, temperature and density are quite different from those of the N 2 - H 2 
case. In this case the velocity and mass fraction profiles are those of a very thin shear layer. 
The temperature and density profiles indicate an only slightly thicker layer with short “tails" 
in the higher speed, colder N 2 stream. Again this behavior is seen in the Lewis and Prandtl 
number profiles of Figures 3 and 6(b). One can see that there is relatively little diffusion of H 2 
from the hotter, slower stream into the N 2 of the colder, faster stream. Note that for this case, 
in contrast to the N 2 - H 2 flow, the density and temperature profiles have only a minimal offset 
from the velocity and mass fraction profiles. The centers of all of the profiles are slightly shifted 
towards the high speed, N 2 side of the layer 

Finally, the velocity, mass fractions, temperature and density profiles, at a free stream Mach 
number of 3, for N 2 — H 2 are shown in Figure 9 and for H 2 - N 2 in Figure 10. The mean profiles 
of velocity, mass fraction and density at this higher Mach number are nearly identical to those 
at M = 0 shown in Figures 7 and 8. This is true for both the case in which N 2 is in the slow 
stream and H 2 is in the fast stream (Figure 9) as well as the reverse (Figure 10). The only 
exception is the temperature profile for the N 2 - H 2 case. The long tail extending into the slow 
speed, N 2 side of the mixing layer is quite different at M — 3 as compared to M = 0. With 
M = 3, the temperature gradient is considerably higher near the center of the layer than at 
M — 0. However, there is almost no difference in the temperature profile between M — 0 and 
M = 3 for the case in which H 2 is in the slow stream and N 2 is in the fast stream. 

5 Conclusions 

By calculating the variation of the Lewis and Prandtl numbers across a mixing layer for all 
thirty binary combinations of H 2 , He, N 2 , 0 2 , Ne and Ar, we have shown that for most pairs of 
gas mixtures there are relatively large variations in the Lewis number across the mixing layer. 
There are also large variations in the Prandtl number for some of these binary gas mixtures. The 
largest variation of both the Lewis and Prandtl numbers occurs for mixtures of the light gases 
He or H 2 with the heavy gases Ar, 0 2 and N 2 . However, the magnitude of the variation is not 
strictly monotone with the molecular weight ratio. Both the Lewis and Prandtl number curves 
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are not symmetric about the center of the mixing layer. This is a result of the non-symmetric 
variation of the temperature across the layer. All of the transport coefficients are functions of 
temperature and so show a corresponding non-symmetric variation across the layer. The “tails 
of the Lewis number and Prandtl number curves are always longest on the side of the mixing 
layer where a lighter gas is diffusing into the higher temperature region as compared to the 
reverse situation. Assuming constant values for the Lewis and Prandtl numbers of a binary gas 
mixture in the shear layer is a poor approximation. 

The general shapes of the velocity, mass fraction, temperature and density profiles for a 
binary mixture of N 2 in the slow stream and H 2 in the fast stream are similar to those obtained 
by Kennedy and Gatski [11] despite their use of an inconsistent definition of the Lewis number. 
This is, in part, because the values of these variables at ± 00 are set by the boundary conditions 
and the equations enforce a smooth transition between these boundary values. The temperature 
and density profiles are shifted towards the N 2 side of the layer as compared to the velocity 
profile. For the opposite case of H 2 in the slow stream and N 2 in the fast stream (not done by 
Kennedy and Gat ski) the profiles of velocity, mass fraction, temperature and density are very 
different from those of the other case. Here the velocity and mass fraction profiles are those of a 
very thin shear layer. The temperature and density profiles are those of a slightly thicker layer 
with quite short “tails” in the higher speed, colder N 2 stream. There is relatively little diffusion 
of H 2 from the hotter, slower stream into the N 2 of the colder, faster stream. For this case, in 
contrast to the N 2 - H 2 flow, the density and temperature profiles have only a minimal offset 
from the velocity and mass fraction profiles. Hence, one would expect that interchanging the 
gases in the slow and high speed streams would have a considerable effect on the flow stability 
and mixing. Finally, for the conditions used in this calculation, the profiles at M = 3 are very 
similar to those for M = 0. 
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Table 1: The ratio of molecular weights W , defined as the molecular weight of the gas at r? = +oc 
divided by the molecular weight of the gas at 77 = - 00 , for the different gases considered in this 
study. The top row corresponds to the gases in the freestream at t? = + 00 , while the first column 
corresponds to the gases in the freestream at 77 = — 00 . 



H 2 

He 

Ne 

n 2 

o 2 

Ar 

H 2 

1.000 

1.986 

10.011 

13.897 

15.873 

19.813 

He 

0.504 

1.000 

5.042 

6.999 

7.994 

9.979 

Ne 

0.100 

0.198 

1.000 

1.388 

1.585 

1.979 

n 2 

0.072 

0.143 

0.720 

1.000 

1.142 

1.426 

o 2 

0.063 

0.125 

0.631 

0.876 

1.000 

1.248 

Ar 

0.050 

0.100 

0.505 

0.701 

0.801 

1.000 
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Lewis Number 


2.6 



Figure 1. Plot of the Lewis number as a function of rj for flu = 0.5, flj = 1.5 and 
Moc = 0. The dashed curves corresponds to the five combinations: He - Ar, H 2 - Ar, 
jVe — Ar, 0 2 — Ar, N 2 — Ar', the solid curves corresponds to the same combinations, 
but interchanged in the free streams. 
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Lewis Number 



11 

Figure 2. Plot of the Lewis number as a function of t] for (3u = 0.5, fir = 1.5 and 
Moo = 0. The dashed curves corresponds to the five combinations: He - 0 2 , H 2 - 02, 
Ne - 02, N 2 - 02, H 2 ~ He', the solid curves corresponds to the same combinations, 
but interchanged in the free streams. 
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Figure 3. Plot of the Lewis number as a function of 7? for 3jj — 0.5, /?r = 1*5 and 
M 0 o = 0. The dashed curves corresponds to the five combinations: He — N 2 , H 2 — A r 2 , 
Ne - N't, He - iVe, H 2 - Ne ; the solid curves corresponds to the same combinations, 
but interchanged in the free streams. 
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Figure 4. Plot of the Prandtl number as a function of 7? for &u = 0.5, /?r = 1.5 and 
A/qo = 0. The dashed curves corresponds to the five combinations: He- Ar , /f 2 - Ar, 
Ne - Ar, 0 2 - -4r, iV 2 - Ar; the solid curves corresponds to the same combinations, 
but interchanged in the free streams. 


16 




Prandtl Number 



Figure 5. Plot of the Prandtl number as a function of rj for flu = 0.5, Pt — 1.5 and 
Mr , o = 0. The dashed curves corresponds to the five combinations: He - 0 % , Hi - Oi , 
Ne - O 2 , Ni - O 2 , Hi - He ; the solid curves corresponds to the same combinations, 
but interchanged in the free streams. 
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Figure 8. Plot of velocity, mass fractions, temperature and density as a function of 
7 / for the binary system H 2 - N 2 . Here, /% = 0.5, fir = 1-5 and M = 0. 
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Figure 9. Plot of velocity, mass fractions, temperature and density as a function of 
r] for the binary system N 2 - # 2 * Here, (5u = 0.5, fa — 1-5 and M <*, = 3. 
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T| 

Figure 10. Plot of velocity, mass fractions, temperature and density as a function 
of Tf for the binary system H 2 - N 2 . Here, 0u = 0.5, f3j = 1.5 and = 3. 
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